APS/123-QED 

Nuclear matter incompressibility coefficient in relativistic and 
nonrelativistic microscopic models 

B. K. Agrawal, S. Shlomo, and V. Kim Au 

Cyclotron Institute, Texas A&M University, College Station, TX 77843, USA 

Abstract 

We systematically analyze the recent claim that nonrelativistic and relativistic mean field (RMF) 
based random phase approximation (RPA) calculations for the centroid energy Eq of the isoscalar 
giant monopole resonance yield for the nuclear matter incompressibility coefficient, Knm, values 
which differ by about 20%. For an appropriate comparison with the RMF based RPA calculations, 
we obtain the parameters of the Skyrme force used in the nonrelativistic model by adopting the same 
procedure as employed in the determination of the NL3 parameter set of the effective Lagrangian 
used in the RMF model. Our investigation suggests that the discrepancy between the values of 
Knm predicted by the relativistic and nonrelativistic models is significantly less than 20%. 

PACS numbers: 21.65+f,24.30.Cz,21.60jz,21.10Re 



1 



I. INTRODUCTION 



The nuclear matter incompressibility coefficient Knm plays an important role in under- 
standing a wide variety of phenomena ranging from giant resonances in finite nuclei to 
heavy-ion collisions and supernova explosions. The measurement of the centroid energy of 
the isoscalar giant monopole resonance (ISGMR) provides a very sensitive method to de- 
termine the value of Knm- Over the last couple of decades, attempts have been made to 
measure very accurately the value of the ISGMR centroid energy, Eq. The recent exper- 
imental data ^ for the Eq in heavy nuclei are accurate enough to provide unambiguous 
information about Knm- However, the theoretical scenario of the Knm still remains unclear. 

In the past, the value of Knm was determined using a macroscopic approach which relies 
on the liquid drop type of expansion for the breathing mode restoring force. It was shown 
that equally good fits can be obtained with values of Knm ranging from 100 — 400 MeV 
In other words, the liquid drop approach can not constrain the value of Knm better than 
50%. On the other hand, the microscopic determination of Knm using the Hartree-Fock 
based random phase approximation (HF-RPA) has undergone a significant improvement 
over time. The earlier nonrelativistic HP— RPA calculations carried out using the Skyrme 
interaction, which nicely reproduced the gross properties of nuclei (such as nuclear binding 
energy and charge radii) and the available data on isovector giant dipole resonance and 
isoscalar giant quadrupole resonance, yielded values of about 370 MeV for the Knm- With 
these interactions, the ISGMR in ^"^Pb was predicted to be located at an excitation energy of 
about 18 MeV. The discovery of ISGMR in ^ospb at an excitation energy of 13.7 MeV ^ led 
to the modification of the Skyrme interaction. Until today, the nonrelativistic calculations 
with Skyrme 0,0] and Gogny interactions predict a value of Knm in the range of 210 — 220 
MeV. We remark here that the long standing problem of the conflicting results deduced for 
Knm from data on the isoscalar giant dipole resonance (ISGDR) and data on the ISGMR 
was explained recently by Shlomo and Sanzhur Q as being due to a missing strength in the 
experimental data for the high energy region of the ISGDR. 

The relativistic mean field based RPA (also referred to as RRPA) calculations, with 
the contribution from negative-energy sea not included, yielded for Knm a value in the 
range of 280 - 350 MeV [gj. Recent RRPA calculations 1^, Q, with the inclusion of 
negative-energy states in the response function, yield a value of Knm between 250 — 270 
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MeV. The discrepancy of about 20% in the value of Knm obtained from relativistic and 
nonrelativistic models is quite significant in view of the accuracy of the experimental data 



currently available on the ISGMR centroid energies. In recent studies 



12, hj it has been 



claimed that these pronounced differences are due to the model dependence of Knm- On 
the other hand, it has been pointed out in Ref. jlj] that the differences in the values of 
Knm obtained in the relativistic and the nonrelativistic models can be attributed, at least in 
part, to the differences in the density dependence of the symmetry energy in these models. 
However, in Ref. , the analysis was restricted to a single nucleus ^osp^ g^^^ ^j^g interaction 
parameters for the several families of the effective Lagrangian considered were fitted only to 
the empirical values of saturation density, binding energy per nucleon in symmetric nuclear 
matter and the charge radius of the ^osp^ nucleus. It may be pointed out here that in Ref. 
jisl ] a reasonable value of Eq for the ^osp]^ nucleus was obtained using an effective force with 
Knm = 400 MeV. But, the same effective force overestimated the value of Eq in case of the 
^°Zr nucleus. Moreover, it has been suggested in Ref. 3| that a wide range of combinations 
of bulk, surface and asymmetry contributions to the finite nucleus compressibility can fit 
the energy of the ISGMR in medium to heavy nuclei. This implies that for a meaningful 
informations about the discrepancy between the relativistic RRPA and the nonrelativistic 
HF-RPA calculations, one must compare the results obtained from these models for several 
nuclei. 

It is noteworthy that in a crude approximation, the uncertainty of about 20% in the 
values of Knm is tantamount to an uncertainty of 10% in the value of Eq. This is because, 
in a semi-classical approach, Eq oc \/Knm- We have shown very recently that the 
calculated value of i?o can deviate by about 5% if the particle-hole space is quite limited 
and/or self-consistency is not properly maintained. We note that in the published literature, 
the calculated values of E^ for ^°^Pb, obtained for the same interaction, differs by up to 0.3 
MeV IHQQ. Therefore, appropriate comparison for the values of Knm obtained from 
the different models is possible only when all the calculations are performed with the same 
procedure and numerical accuracy. 

In this work we take a close look at the issue of the model dependence of the nuclear 
matter incompressibility coefficient derived from the ISGMR centroid energy. Toward this 
purpose, we generate different parameter sets for the Skyrme interaction and perform highly 
accurate calculations for the ISGMR strength function for several nuclei using the HF based 
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RPA approach. For the sake of true comparison, the calculations using different parameter 
sets of Skyrme interaction are performed following exactly the same numerical procedure. 
The values of the Skyrme parameters are obtained by a least square fit to exactly the same 
experimental data for the nuclear binding energies, charge radii and neutron radii as adopted 
in Ref. in determining the NL3 parameter set for an effective Lagrangian used in the 
RMF model. We find that the model dependence is rather weak and the differences in the 
values of Knm in the relativistic and the nonrelativistic models essentially arise from the 
differences in the nuclear matter properties, in particular, in the values of the symmetry 
energy coefficient (J), associated with these models. 



II. FORMALISM 



In a self-consistent HF-RPA calculation [20|, one starts by adopting a specific effective 



nucleon-nuc^ 
the form 



eon interaction, V12. In this work we shall use the Skyrme type interaction of 
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to (1 + XoPfs) 5(ri - rs) + -h (1 + x^P'^^) 



X 



k \^6{vi - V2) + 5(ri - ra) kl^ + (1 + X2P12) k i2<^(ri - rs) k 

r2 
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+ it3(l + X3Pf2)p"(^ 



5{vi - V2) 



-\-iWo k i2(^(ri - r2)(at + 0^) x k 12 



(1) 



where, is the spin exchange operator, is the Pauli spin operator, k 12 = — z( Vi — V2)/2 
and k 12 = — z(Vi — V2)/2. Here, right and left arrows indicate that the momentum 
operators act on the right and on the left , respectively. The parameters of the Skyrme force 
are obtained by fitting the HF results to the experimental data for the bulk properties of 
finite nuclei. Once the HF equations are solved using an appropriate parameter set for the 
Skyrme interaction, then one obtains the RPA Green's function 



G — Go{l + Vp^hGo 



(2) 



where, Vp-h is the particle-hole {p — h) interaction consistent with V12 and Gq is the free 
p — h Green's function . For the single-particle operator 



(3) 



i=l 



the strength function is given by 

SiE) = y2\{0\F\n)\'6iE- E^) = -Im [Tr (fGf)] . (4) 

n 

The steps involved in the relativistic mean field based RPA calculations are analogous 
to those described above for the nonrelativistic HF-RPA approach. However, the nucleon- 
nucleon interaction in case of relativistic mean field models are generated through the ex- 
change of various mesons. An effective Lagrangian which represent a system of interacting 
nucleons looks like [l^, 

1 2 

C = i){'-i{id- g^u - Qpp T ~eA) - m- g^a) i' + ^ (da) 

+ \ml-p' - \f,.F^' (5) 

which contains nucleons with mass m; a, uj, p mesons; the electromagnetic fields; and 
non-linear self-interactions of the a field, 

U{cr) = \mla' + ^-g,a' + \g,a\ (6) 

The Lagrangian parameters are usually obtained, as in the case of nonrelativistic mean field 
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calculations, by fitting procedure to some bulk properties of a set of spherical nuclei 
The values of various coupling constants and the meson masses appearing in Eqs. (0) and 
© for the most widely used parameter set NL3 are mo- = 508.194 MeV, = 782.501 
MeV, nip = 763.000 MeV, g^ = 10.217, g^ = 12.868, c/p = 4.474, g2 = -10.431 fm'^ and 
^3 = -28.885. 

III. RESULTS AND DISCUSSIONS 
A. HF results 

In the present work, for an appropriate comparison, we carry out a least square fit to 
the same experimental data for the nuclear binding energies, charge radii and neutron radii 
used in Ref. jl^ to obtain the NL3 parameter set. Furthermore, we deal with the centre of 
mass correction to the total binding energy, finite size effects of the proton and the Coulomb 
energy in the way similar to that employed in determining the NL3 parameter set in Ref. 



[igj l . It may be pointed out that pairing is not included in our HF calculations and instead of 
the 10 nuclei considered in Ref. 2, we consider 7 closed shell nuclei. The open shell nuclei 
^^Ni, ^^"^Sn and ^^^Pb are excluded from our least square fit. We also ignore the proton 
and the neutron pairing gaps in case of ^°Zr and ^^^Sn nuclei, respectively. However, we 
have verified that if we increase the error bars for the experimental data on these two nuclei 
in order to compensate for the pairing, we find that the values of the Skyrme interaction 
parameters remains practically unaltered. 

Since, the main objective of this paper is to delineate the differences in the value of 
Knm predicted by the relativistic and nonrelativistic mean field based RPA calculations, we 
generate a Skyrme interaction having Knm and J very close to those associated with the 
NL3 parameter set, i.e., 271.76 and 37.4 MeV, respectively. Also, as most of the calculation 
claiming the model dependence are restricted to the ISGMR centroid energy for the single 
heavy nucleus ^°^Pb, we generate a parameter set by demanding a very high accuracy for 
the root mean square charge radius of ^°^Pb. We denote this parameter set as SK272. In 
addition, we also obtain a parameter set SK255 having characteristics very much similar to 
the SK272 parameter set, but, Knm is taken to be 255 MeV. 

In Table U we give the values of the Skyrme parameters SGII together with the new 
parameter sets SK272 and SK255. In Table ITTl we compare the nuclear matter properties for 
the SK272 and SK255 interactions with the corresponding ones obtained from the NL3 and 
SGII interactions. The quantities po? m* /m and L in this table denote the saturation density, 
effective nucleon mass and the slope of the symmetry energy coefficient (L = SpodJ/dpo), 
respectively. In column 3 of Table IIIII we have given the experimental data for the total 
binding energy, E, charge radii, r^, and neutron radii, r„, for the nuclei used in Ref. [3| 
for determining the NL3 parameter set and in our fit, together with the assumed error bars 
in percent. The values obtained from the parameter sets SK272 and SK255 are shown in 
columns 5 and 6, respectively. For the sake of comparison, in Table IIIII we also give in 
columns 4 and 7 the results for the NL3 and SGII interactions, respectively. It is evident 
from this table that the quality of our fit to the experimental data is quite comparable to 
the results obtained with the NL3 and SGII interactions. We have also listed the values of 
Ar = r„ — Tp, the difference between rms radii for neutrons and protons, (not included in 
the fit). Experimental values for Ar are obtained from the data for the neutron and charge 



radii and by using = \/r1 — 0.64. It is interesting to note that the values of Ar obtained 
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for the SK272, SK255 and NL3 interactions are closer and are quite large compared with 
the corresponding results for the SGII interaction. This is consistent with the fact that the 
values of the slope of the symmetry energy (L) associated with the SK272, SK255 and NL3 
interactions are significantly larger than that associated with the SGll interaction (see Table 
inj- Since, as is well known j2^, the value of Ar is sensitive to the density dependent form 
adopted for the symmetry interaction. 



B. RPA results 

We have demonstrated very recently that the strength function for giant resonances 
are quite sensitive to the various numerical approximations made. By numerical approxima- 
tions we essentially mean, the size of the box used for the discretization of the continuum, 
restriction imposed on the maximum energy for the particle- hole excitations {E'^^) , smear- 
ing parameter (r/2) used to smear the strength function, etc. We have shown in Ref. 
that in order to reproduce the results obtained in the continuum RPA calculation the size 
of the box should be consistent with the value used for the smearing width. For r/2 = 1 
MeV, one must use a large box of size 72 fm. In this work, we have used a box of size 90 
fm and 1/2 = 1 MeV. From Ref. Q we also note that in order to obtain a reliable value of 
the ISGMR centroid energy, Eq, accurate within 0.1 MeV, one must use E^^ > 400 MeV. 
Here, the centroid energy is given by Eq = mi /mo, where mo and mi are the non-energy- 
weighted and energy- weighted sums of S{E) of Eq. (jU, respectively. In the present work, 
the lowest value of E^^ we have used is higher than 500 MeV. In addition, the centroid 
energy depends strongly on the range adopted for the excitation energy interval of the giant 
resonance. Nevertheless, one often encounters in the published literature that the values of 
the centroid energies are given without any reference to the corresponding excitation energy 
range considered. For instance, in the case of the ^osp]^ nucleus we find for the SGII pa- 
rameter set that Eq = 13.7, 13.9, 14.4 MeV for the excitation energy ranges — 40, — 60 
and 10 — 40 MeV, respectively. These differences are quite significant, since, as pointed out 
earlier, a variation of 5% in Eq corresponds to a change in Knm by 10%. In what follows, we 
shall concentrate mainly on the results for Eq obtained by integrating the strength function 
over the energy range — 60 MeV, since the RMF results presented in Ref. for the 
NL3 parameter set were obtained using the same energy range and the strength function 
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was smeared using r/2 = 1 MeV 



2^. 



In Table HVl we give the results for the ISGMR centroid energy obtained using the param- 
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for the 



eter sets SK272 and SK255 and compare them with the RRPA results of Ref. 
NL3 interaction. It is clear from this table that the parameter set SK272 yields values for 
Eq which are higher by about 4 — 6% than those obtained using the NL3 parameter set in 
the RRPA calculations. This implies that if we reduce the compressibility by about 10%, we 
can reproduce reasonably well the RRPA results for the NL3 interaction. For this reason we 
have generated another parameter set SK255 with Knm = 255 MeV, keeping J = 37.4 MeV 
(see Table HH). In fact, we see that the differences between the values of the Eq obtained 
from the parameter sets SK255 and NL3 are on the level of the uncertainty associated with 
the experimental data for Eq. We emphasize that though the values of Knm associated with 
the SGII and SK255 parameter sets differ by about 40 MeV, the values of Eq for ^osp]^ 
nucleus for these interactions are close within 0.3 MeV. Thus, by requiring a fixed value of 
Eq we find that an increase in J by 10% leads to an increase in Knm by about 5%. We also 
compute the values of Enover the same energy range as used in experimental determination 
of the centroid energy jl|. It should be noted that the experimental values of the energy 
range, given in the Table HVi are more or less symmetric around the corresponding Eq. It 
can be seen from Table HV| that for the parameter set SK255, we obtain a good agreement 
with the experimental data for Eq, calculated over the experimental excitation energy range. 
We remark that, in our calculations for the ^°^Pb nucleus with the SK272 and SK255 pa- 
rameter sets, the peak energy for the isovector giant dipole resonance is 13.2 and 13.3 MeV 
respectively, which is in good agreement with the experimental value of 13.3 ±0.1 MeV 

C. Conclusions 

In summary, we have analyzed in detail the recent claim that the nuclear matter incom- 
pressibility coefficient K^m extracted from the ISGMR centroid energy calculated using the 
relativistic and nonrelativistic based RPA models differ by about 20%. For a meaningful 
comparison, we have generated parameter sets for the Skyrme interaction by a least square 
fitting procedure u sing exactly the same experimental data for the bulk properties of nuclei 
considered in Ref. 19| for determining the NL3 parameterization of an effective Lagrangian 
used in the relativistic mean field models. Further, we also demanded in our fitting proce- 
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dure that the values of Knm, J and the charge radius of the ^°^Pb nucleus should be very 
close to the results obtained with the NL3 interaction. The parameter sets thus obtained 
were used to calculate the ISGMR centroid energy for several nuclei. For the parameter set 
SK272 {Knm = 272 MeV), the calculated values of Eq are higher by about 5% compared to 
the corresponding NL3 results. This implies that the difference in the value of Knm obtained 
in the relativistic and the nonrelativistic microscopic models could be at most 10%. In view 
of this, we generate another parameter set having Knm — 255 MeV. As expected, the pa- 
rameter set associated with K^m — 255 MeV, yields for the ISGMR centroid energies values 
which are quite close to the NL3 results. Moreover, for the SK255 parameter set, we find a 
good agreement with experimental data for Eq for all the nuclei considered, provided, the 
corresponding excitation energy ranges used in determining Eq arc the same as those used 
in obtaining the experimental data. We have thus shown that the difference in the values of 
Knm obtained in the relativistic and nonrelativistic models is mainly due to the differences 
in the values of the symmetry energy coefficient (J) and its slope (L) associated with these 
models. 



[1] D. H. Youngblood, H. L. Clark and Y. W. Lui, Phys. Rev. Lett. 82, 691 (1999). 

[2] J. M. Pearson, Phys. Lett. B271, 12 (1991). 

[3] S. Shlomo and D. H. Youngblood, Phys. Rev. C 47, 529, (1993). 

[4] N. Marty, et al, Nucl. Phys. A230, 93 (1975). 

[5] G. Colo, P. F. Bortignon, N. Van Giai, A. Bracco and R. A. Broglia, Phys. Lett. B276, 279 
(1992). 

[6] L Hamamoto, H. Sagawa and H. Z. Zhang, Phys. Rev. C 56, 3121 (1997). 

[7] J. P. Blaizot, J. F. Berger, J. Decharge and N. Girod, Nucl. Phys. A591, 435 (1995). 

[8] S. Shlomo and A. I. Sanzhur, Phys. Rev. C 65, 044310 (2002); S. Shlomo, Pramana-J. 57, 
557 (2001). 

[9] Zhongyu Ma, Nguyen Van Giai and Hiroshi Toki, Phys. Rev. C 55, 2385 (1997). 

[10] Zhong-yu Ma, Nguyen Van Giai, A. Wandelt, D. Vretenar, Nucl. Phys. A686, 173 (2001). 

[11] D. Vretenar, T. Niksic and P. Ring, Phys. Rev. C, in press. 

[12] Nguyen Van Giai, P. F. Bortignon, G. Colo, Zhongyu Ma and M. Quagha, Nucl. Phys. A687, 



9 



44c (2001). 

[13] T. Niksic, D. Vretenar and P. Ring, Phys. Rev. C 66, 064302 (2002). 
[14] J. Piekarewicz, Phys. Rev. C 66, 034305 (2002). 

[15] F. E. Serr, Proc. Nuclear Physics Workshop, Trieste (Italy), October 1981, eds. C. H. Dasso, 

R. A. Broglia and A. Winthelr (North- Holland, Amsterdam, 1982). 
[16] B. K. Agrawal, S. Shiomo and A. I. Sanzhur, Phys. Rev. C 67, 034114 (2003). 
[17] Gianluca Colo, Nguyen Van Giai Pier F. Bortignon and Maria R. Quaglia, RIKEN Review 

No. 23 (July 1999). 

[18] Gianluca Colo, Nguyen Van Ciai Pier F. Bortignon and Maria R. Quaglia, Physics of Atomic 

Nuclei Vol. 64, No. 6, 1044 (2001). 
[19] G. A. Lalazissis, J. Konig and P. Ring, Phys. Rev. C 55, 540 (1997). 
[20] G. F. Bertsch and S. F. Tsai, Phys. Rep. 18, 125 (1975). 
[21] D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626, (1972). 
[22] P. -G. Reinhard, Rep. Prog. Phys. 52, 439 (1989). 
[23] E. Friedman and S. Shiomo, Z. Physik A283, 67 (1977). 
[24] Nguyen Van Giai and H. Sagawa, Phys. Lett. B106, 379 (1981). 
[25] Zhong-yu Ma, private communication. 
[26] J. Ritman et al, Phys. Rev. Lett. 70, 533 (1993). 



10 



TABLE I: Skyrme parameters for different interactions used in the present calculations. Value of 
the parameters for the SGII interaction are taken from Ref. [2^. 



Parameter 


SK272 


SK255 


SGII 


io(MeV • fm^) 


-1496.84 


-1689.35 


-2645 


ii(MeV • fm^) 


397.66 


389.30 


340 


t2(MeV • fm^) 


-112.82 


-126.07 


-41.9 


t3(MeV • fm^(^+")) 


10191.64 


10989.60 


15595 


Xn 


0008 


-0.1461 


0.09 


xi 


0.0102 


0.1160 


-0.0588 




0.0020 


0.0012 


1.425 




-0.5519 


-0.7449 


0.06044 


alpha 


0.4492 


0.3563 


1/6 


VFo(MeV • fm^) 


106.58 


95.39 


105 



TABLE II: Nuclear matter properties calculated from RMF theory with the NL3 parameter set 
and the nonrelativistic HF calculations with different Skyrme parameter sets. The "experimental 
data" are the ones used in Ref. jl9| in the least square fit together with the bulk properties for 
finite nuclei in obtaining the NL3 parameter set. The values in parenthesis represent the error bars 
(in percent) used in the fit. 





Exp. 


NL3 


SK272 


SK255 


SGII 


E/A{MeY) 


-16.0(5) 


-16.299 


-16.280 


-16.334 


-15.67 


Knm (MeV) 250.0(10) 


271.76 


271.55 


254.96 


214.57 


Po (fm 3) 


0.153(10) 


0.148 


0.155 


0.157 


0.159 


m* /m 




0.60 


0.77 


0.80 


0.79 


J (MeV) 


33.0(10) 


37.4 


37.4 


37.4 


26.8 


L (MeV) 




118.5 


91.7 


95.0 


37.6 
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TABLE III: Experimental data for the total binding energy (E ) in MeV, charge (rc) and neutron 
(r„) radii in fm for the nuclei used in the least square fitting procedure. The parameter sets thus 
obtained are named as SK272 and SK255. For comparison we also give the results obtained from 
the NL3 and SGII interactions. The values in parenthesis represent the error bars (in percent) 
used in the fit. The quantity Ar = r„ — rp (not included in the fit) is the difference between the 
radii for the neutrons and protons. 



Nucleus Property 


Exp. 


NL3 


SK272 


SK255 


SGII 


160 


E 


-127.62(0.1) 


-128.83 


-127.76 


-128.05 


-131.93 




' c 


2.730(0.2) 


2.730 


2.800 


2.813 


2.793 








2.580 


2.662 


2.674 


2.650 


40Ca 


E 


-342.06(0.1) 


-342.02 


-341.35 


-342.50 


-342.42 




' c 


3.450(0.2) 


3.469 


3.496 


3.504 


3.490 




' 71 


3.370(2.0) 


3.328 


3.363 


3.369 


3.348 




Ar 


0.014 


-0.047 


-0.041 


-0.043 


-0.049 


48Ca 


E 


-416.00(0.1) 


-415.15 


-414.17 


-413.89 


-418.22 




' c 


3.451(0.2) 


3.470 


3.524 


3.531 


3.526 




' n 


3.625(2.0) 


3.603 


3.635 


3.649 


3.582 




Ar 


0.268 


0.227 


0.203 


0.210 


0.147 


9°Zr 


E 


-783.90(0.1) 


-782.63 


-782.73 


-783.28 


-775.49 




' c 


4.258(0.2) 


4.287 


4.282 


4.286 


4.286 




' n 


4.289(2.0) 


4.306 


4.310 


4.317 


4.266 




Ar 


0.107 


0.094 


0.103 


0.106 


0.056 




E 


-988.69(0.1) 


-987.67 


-982.37 


-984.48 


-971.66 




rc 


4.627(0.2) 


4.611 


4.617 


4.619 


4.630 






4.692(2.0) 


4.735 


4.696 


4.701 


4.639 




Ar 


0.135 


0.194 


0.149 


0.152 


0.079 


Sn 


E 


-1102.90(0.1) 


-1105.44 


-1097.36 


1100.04 


-1105.17 




rc 




4.709 


4.725 


4.726 


4.735 








4.985 


4.964 


4.975 


4.867 


208pj-, 


E 


-1636.47(0.1) 


-1639.54 


-1631.78 


-1637.48 


-1622.21 




Tc 


5.503(0.2) 


5.520 


5.503 


5.503 


5.519 






5.593(2.0) 


125.741 


5.687 


5.694 


5.597 




Ar 


0.118 


0.279 


0.213 


0.250 


0.13(3 



TABLE IV: ISGMR centroid energy Eq = mi/niQ (in MeV) obtained by integrating over the energy 
range uJi — L02 MeV with the strength function smeared by using r/2 = 1 MeV. The experimental 
data is taken from Ref. \m. 



Nucleus 


^1 




Exp. 


NL3 SK272 SK255 SGII 


90 Zr 


- 


- 60 




18.7 


20.0 


18.9 


18.3 




10 


- 26 


17.89±0.20 




19.3 


18.4 


17.9 




- 


-60 




17.1 


18.0 


17.5 


16.6 




10 


- 23 


16.07±0.12 




17.4 


16.9 


16.3 




- 


-60 




16.1 


17.1 


16.4 


15.6 




10 


- 22 


15.39±0.28 




16.5 


15.9 


15.2 


208p]g 


- 


-60 




14.2 


14.7 


14.2 


13.9 




8 - 


-21 


14.17±0.28 




14.2 


13.8 


13.6 
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